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We investigate through direct molecular mechanics calculations the geometrical properties of 
hydrocarbon mantles subjected to percolation disorder. We show that the structures of mantles 
generated at the critical percolation point have a fractal dimension df « 2.5. In addition, the 
solvent access surface As and volume Vs of these molecules follow power-law behavior, As ~ L"^ and 
Vs ~ L"^ , where L is the system size, and with both critical exponents aA and ay being significantly 
dependent on the radius of the accessing probing molecule, rp. Our results from extensive simulations 
with two distinct microscopic topologies (i.e., square and honeycomb) indicate the consistency of the 
statistical analysis and confirm the self-similar characteristic of the percolating hydrocarbons. Due 
to their highly branched topology, some of the potential applications for this new class of disordered 
molecules include drug delivery, catalysis, and supramolecular structures. 
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In recent years there has been a great scientific and 
technological advance in the design and development of 
new materials. The use of molecular modelling tech- 
niques is certainly responsible for the significant contri- 
bution to this activities of most theoretical and compu- 
tational studies Ij. New improved catalysts, original re- 
active processing and the chemistry and physics of poly- 
mers are just a few among several examples representing 
the multitude of research fields that have been success- 
fully approached employing molecular modelling tools. 
For instance, an important research work with poten- 
tial applications in the field of polymer science refers to 
the real possibility of manipulating and controlling the 
properties of complex macromolecules through adequate 
polymerization routes and specific chemical processing 
strategies. 

An inspiring example where molecular modelling plays 
an important role in polymer science is the design 
and characterization of functional polymers. By defini- 
tion, the term functional refers to the property of some 
molecules that contain reactive end-groups to preserve 
the integrity of their molecular backbone while partic- 
ipating in chemical reactions. Because dendrimers are 
highly branched macromolecules and, therefore, have a 
large number of end-groups, they represent a natural 
class of candidates for substrates of functional polymers. 
Among many possible applications for dendrimers, we 
mention their potential use for drug delivery 0| , and the 
ability to act as catalysts supports [j] and supramolecular 
structures In a recent theoretical study, the molecu- 
lar topology and surface accessibility of dendrimers have 
been investigated through extensive molecular mechanics 
and molecular dynamics simulations 5]. Accordingly, it 
has been shown that all dendrimcr families exhibit non- 
spherical and shape-persistent structures, with a molec- 
ular fractal dimension close to 2.5 and an invariant sur- 



face fractal dimension approximately equal to 2, regard- 
less the number of generations (branches) composing the 
molecular geometry. 

In the present letter, we propose a new class of macro- 
molecules that are also highly branched, but possess in- 
trinsic disorder in their compositional and structural con- 
formation. Here we use percolation as a model for dis- 
order. Due to its simplicity, the percolation theory 0,0 
has been extensively applied to represent disordered ma- 
terials in several research fields of scientific and tech- 
nological relevance, such as condensed matter physics 
8], flow through porous media 9], and heterogeneous 
chemistry [wl 111 III jT| , among others. More precisely, 
percolation is a purely geometric model in the sense that 
we have to populate all lattice sites (or bonds) with a 
prescribed probability p. For small values of p, only fi- 
nite clusters are present. By increasing p, one can find a 
threshold or a critical value Pc for which an "infinitely" 
connected object, the spanning cluster, is generated. In 
fact, the critical geometry of the spanning cluster [3. ITU 
is an example of random fractal that has been frequently 
used as a representation for real disordered systems. The 
aim here is to investigate through molecular mechanics 
simulations the geometrical properties of hydrocarbon 
molecules generated from carbon mantles subjected to 
site percolation disorder. 

Starting from a square or an hexagonal lattice of size L, 
we extract the spanning cluster and assign to each site a 
carbon atom that is connected by single bonds to its near- 
est neighbors. The valence of each carbon in the cluster 
is then adjusted to 4 by adding the necessary hydrogen 
atoms, so that, in average, a hydrocarbon with fractional 
stoichiomctry is produced, namely, C^^Hj,. The bond 
length in the initial configuration of this molecule is set 
to an arbitrary guess value (e.g., lA) and a molecular 
mechanics method is applied to optimize the molecular 
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geometry 'l^. Molecular mechanical force fields use the 
equations of classical mechanics to describe the potential 
energy surfaces and physical properties of molecules. Un- 
der this framework, a molecule is viewed as a collection 
of atoms that interact with each other by simple analyt- 
ical functions. Here we apply the standard MM+ force 
field [31 to obtain the optimized geometry for each real- 
ization of the percolating hydrocarbon. For comparison, 
a limited number of computational simulations have also 
been performed with a more realistic approach, namely, 
the PM3 semi-empirical and self-consistent-field molec- 
ular orbital method For all practical purposes, the 
relative differences found between the geometrical prop- 
erties calculated here resulting from the PM3 technique 
and those obtained using the MM+ molecular mechan- 
ics are always smaller than 3%. 

In Fig. la we show the the original square lattice 
of carbon sites (planar) and the optimized geometries 
computed for a carbon mantle that is fully occupied 
(p = 1). When submitted to the optimization proce- 
dure, the molecule bends spontaneously to assume the 
shape of a corrugated tile. In the case of an hexagonal 
lattice, the resulting molecular structure remains planar. 
The structures shown in Figs, lb and Ic correspond to 
typical configurations of hydrocarbons mantles generated 
from square and hexagonal lattices, respectively, but now 
with the probability p set right above the critical perco- 
lation point Pc- 

We performed simulations for square and hexagonal 
lattices generated at the critical point {pc ~ 0.593 and 
0.696, respectively), and with linear sizes ranging from 
L = 10 to 32. In the case of the square lattice, for each 
linear size L = 12, 15, 18, 21, 24, 27, and 30, we used 300, 
300, 150, 150, 80, 80, and 50 realizations, respectively. 
For the hexagonal geometry, the calculations have been 
carried out with 300, 150, 150, 80, 80, and 50 reahza- 
tions for L = 10, 14, 18, 22, 26, and 30, respectively For 
each realization, once the convergence for the minimum 
energy structure is numerically attained, we compute a 
set of relevant geometrical properties that are then aver- 
aged over all samples of same size L. For instance, the 
radius of gyration Rg represents a meaningful measure 
to characterize objects of complicated geometry 0]. It 
is defined here as i?^ = Yliti ~ ^of/Na, where Na is 
the number of atoms, ro = SilTi ri is the position 

of each atom i in the molecule, regardless of its type, i.e., 
hydrogen or carbon. The results shown in Fig. 2 clearly 
indicate that the radius of gyration displays a typical 
power-law behavior for both networks topologies, 

Rg - M^/'^f , (1) 

where M is the average molecular weight and d/ is the 
fractal dimension of the disordered molecules. The best 
fits to the data yield exponents that are statistically sim- 
ilar, namely, 0.40 ± 0.03 and 0.38 ± 0.02 corresponding 
to the fractal dimensions df — 2.50 and 2.63, for the 
square and hexagonal lattices, respectively. Considering 
the small sizes of the lattices employed in the simulations 




FIG. 1: (a) Hydrocarbon mantle of size L = 18 generated 
from a fully occupied square lattice {p — 1). The hydrocarbon 
chains on the left are the initial two-dimensional structures 
used in the optimization procedure, while those on the right 
are the corresponding three-dimensional topologies resulting 
from the molecular mechanics simulations, (b) The same as 
in (a), but for a value of p right above the critical percolation 
point Pc- (c) The same as in (b), but for an hexagonal lattice 
topology. 



due to computational limitations, the agreement between 
these two exponents suggests that the complex geometry 
of percolating hydrocarbons may constitute a single class 
of universality. 

At the microscopic scale it is important to determine 
the volume and surface of a molecule. Here we use the 
concept of van der Waals volume (surface) which is the 
sum of the volume (surface) of each atom composing a 
molecule, calculated as a sphere with the corresponding 
van der Waals radius. By definition, the van der Waals 
radii are simply the radii of the spheres connecting two 
non-bonded atoms Q . The inset of Fig. 2 shows the log- 
log plot of the van der Walls volume Vw against the size 
L of the system for both square and hexagonal lattices. 
From the least-square fits to the two data sets we observe 
that the results for both lattice topologies can be well 
described by the scaling relation, Vw ^ L'^^ , where dw 
is the critical exponent. If we consider that each carbon 
gives approximately the same contribution to the overall 



3 



1.30 




surface 



FIG. 2: Logarithmic plot showing the dependence of the ra- 
dius of gyration Rg on the average molecular weight M of 
percolating hydrocarbons generated from square (circles) and 
hexagonal (triangles) lattices. The straight lines are the least- 
square fits to the data with the numbers indicating the slopes, 
0.40 ± 0.03 (circles) and 0.38 ± 0.02 (triangles). The inset 
shows the log-log plot of the van der Waals volume against 
the linear dimension of the system for square (circles) and 
hexagonal (triangles) lattices. 



volume of the molecule, the values dw — 1-84 ± 0.03 
and 1.87 ± 0.02 obtained for the square and hexagonal 
lattice, respectively, can be readily compared with the 
fractal dimension of the spanning cluster at the critical 
percolation point, which is ~ 1.89 in two dimensions 0. 
The small discrepancies found here are due to finite-size 
effects as well as the presence of the hydrogens atoms 
necessary to balance the valence of the carbon sites. 

As already mentioned, the investigation of molecular 
volume and surface is essential for the development of 
many chemical processes of scientific and technological 
relevance. Indeed, the molecular architecture can play a 
major role in determining the potential application and 
efficiency of functional macromolecules jj, jig • In rela- 
tion to that, the concepts of solvent accessible surface Ag 
and volume Vs jT3| have been extensively used to charac- 
terize and compare the geometrical aspects of the interac- 
tion between a molecule and different types of solvent. As 
depicted in Fig. 3, the Ag is obtained by rolling a spher- 
ical probe of diameter 2rp corresponding to the size of 
the solvent (e.g., water or ethanol) on the van der Waals 
surface (Aw) of the molecule. The Vg is the volume for 
which the boundary is the Ag. Since macromolecules like 
the percolating hydrocarbons studied here contain small 
gaps, pockets and clefts which are sometimes too small to 
be penetrated even by a solvent molecule like water, the 
Ag becomes gradually smoother as the size of the solvent 
used for probing increases. In the limiting case where 
is set to zero, we recover the Aw of the molecule. 

In Fig. 4 we show the log-log plot of the average Vg 
against the linear size of percolating hydrocarbons gen- 



FIG. 3: The solvent accessible surface and volume of each re- 
alization of a percolating hydrocarbon are obtained by rolling 
a spherical probe of radius Vp on its van der Waals surface. 
For each lattice topology (square or hexagonal) and a given 
value of the linear dimension L, the average surface As and 
volume Vs are calculated from several realizations of the dis- 
ordered molecules. 
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FIG. 4: Log-log plot of the variation of the average accessible 
volume Vs with the linear dimension L for four different values 
of the solvent probe radius rp. The inset shows the same vari- 
ation, but for the average accessible surface As- The straight 
lines are the best fits to the data sets of the scaling relations 
As oc L"^ and Vs oc L"^ , with the slopes corresponding to 
the exponents aA and av ■ These results are for percolating 
hydrocarbons generated from a square lattice topology. 



crated from square lattices, and calculated for four dif- 
ferent values of rp. These results provide clear evidence 
that the Vg of such disordered molecules increases as a 
power-law, 
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FIG. 5: Dependence of the scaling exponents aA and av 
on the solvent probe radius Vp for percolating hydrocarbons 
extracted from a square lattice. 



with an exponent ay that significantly depends on the 
size Tp of the probing molecule. The inset of Fig. 4 shows 
that the Ag also follows a power-law, Ag ^ L"^ , with an 
exponent aA that is, like in the Vs case, nonuniversal with 
respect to the solvent radius. We find an entirely similar 
behavior for both Vs and As obtained from molecular 
mechanics simulations performed with critical hexagonal 
lattices. 

Figure 5 shows the dependence of ay and ua on the ra- 
dius Tp for the case of the square topology. Starting from 
values that are close to the fractal dimension of the span- 
ning cluster in two-dimensions, both exponents decrease 
monotonically with r^, reflecting the gradual smoothness 
of the As as the size of the solvent molecule increases. In 
all cases, we find a maximum relative deviation of 5% 
between scaling exponents of square and hexagonal lat- 



tices computed at the same value of r^. As expected, 
these systems belong to the same universality class, de- 
spite the difference in the details of their microscopic ge- 
ometry as well as the fact that a larger contribution of 
hydrogen atoms is needed to balance the critical hexag- 
onal molecules, in comparison to their square counter- 
parts of same size. These results indicate that any type 
of physical or chemical interaction between a percolating 
hydrocarbon and a given molecule can be dramatically 
dependent on the molecule size. For instance, even if 
uniformly functionalized, these highly structured hydro- 
carbons may display enormous discrepancies in reactivity 
if the reagents have slightly different sizes. 

In conclusion we have studied the average geometrical 
features of molecules generated from hydrocarbon man- 
tles subjected to percolation disorder. Our results show 
that these nanostructures are self-similar and display a 
rich variety of scaling behavior in their solvent accessi- 
ble volume and surface properties. We expect these dis- 
ordered molecules to be experimentally feasible. As an 
important step in this direction, ultrathin sheet-like car- 
bon nanostructures, also called carbon nanosheets, have 
been recently synthesized under high controllability by 
a plasma enhanced process of chemical vapor deposition 
|16|. Among other characteristics, these molecules have a 
very high surface-to- volume ratio and sharp edges which 
are attractive for fuel cells and micro-electronic technolo- 
gies. Under a different perspective, Andrade et al. 
proposed a growth mechanism for branched polymers 
where self-organization leads the system spontaneously 
to a percolation-like critical state that is entirely simi- 
lar to the disordered molecules investigated here. In this 
study, a conceivable experimental scheme is suggested for 
the implementation of the method that could lead to the 
generation of percolation hydrocarbons. 

We thank the Brazilian agencies CNPq, CAPES and 
FUNCAP for financial support. 



[1] A. R. Leach, Molecular Modelling: Pnnctples and AppU- 

cations (Pearson, Edinburgh, 2001). 
[2] R. Esfand and D. A. Tomalia, Drug Discovery Today 6, 

427 (2001). 

[3] H. Brunner, J. Organomet. Chem. 500, 39 (1995). 

[4] T. Emrick and J. M. J. Frechet, Curr. Opin. Coll. Inter- 
face Sci. 4, 15 (1999). 

[5] S. Pricl, M. Fermeglia, M. Ferrone, and A. Asquini, Car- 
bon 41, 2269 (2003). 

[6] D. Stauffer and A. Aharony, Introduction to Percolation 
Theory (Taylor Francis, Philadelphia, 1994). 

[7] M. Sahimi, Applications of Percolation Theory (Taylor & 
Francis, London, 1994). 

[8] T. Nakayama and K. Yakubo, Fractal Concepts in Con- 
densed Matter Physics (Springer- Verlag, Berlin, 2003). 

[9] M. Sahimi, Flow and Transport in Porous Media and 
Fractured Rock 
(VCH, Boston, 1995). 
[10] D. Avnir (Ed.) (1989) The Fractal Approach to Hetero- 



geneous Chemistry: Surfaces, Colloids, Polymers (Wiley, 
Chichester, 1989). 

[11] S. Havlin and A. Bunde, "Percolation II: Transport in 
Percolation Clusters" in Fractals and Disordered Sys- 
tems, Second Edition edited by A. Bunde and S. Havlin 
(Springer- Verlag, Berlin, 1996). 

[12] J. S. Andrade Jr., D. Street, Y. Shibusa, S. Havlin, and 
H. E. Stanley, Phys Rev. E 55, 772 (1997). 

[13] M. H. A. S. Costa, A. D. Araujo, H. F. da Silva, and J. 
S. Andrade Jr., Phys . Rev. E 67, 061406 (2003). 

[14] N. L. AUinger, J. Am. Chem. Soc. 99, 8127 (1977). 

[15] J. J. P. Stewart, J. Comput. Chem. 10, 209 (1989). 

[16] J. Wang, M. Zhu, R. A. Outlaw, X. Zhao, D. M. Manos, 
and B. C. HoUoway, Carbon 42, 2867 (2004). 

[17] M. L. Connolly, J. Am. Chem. Soc. 107, 1118 (1985). 

[18] J. S. Andrade Jr., L. S. Lucena, A. M. Alencar, J. E. 
Freitas, Physica A 238, 163 (1997). 



